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ABSTRACT 

Distances to stars are key to revealing a three-dimensional view of the Milky Way, yet their 
determination is a major challenge in astronomy. Whilst the brightest nearby stars benefit 
from direct parallax measurements, fainter stars are subject of indirect determinations with 
uncertainties exceeding 30%. We present an alternative approach to measuring distances using 
spectroscopically-identified twin stars. Given a star with known parallax, the distance to its 
twin is assumed to be directly related to the difference in their apparent magnitudes. We 
found 175 twin pairs from the ESO public HARPS archives and report excellent agreement 
with Hipparcos parallaxes within 7.5%. Most importantly, the accuracy of our results does not 
degrade with increasing stellar distance. With the ongoing collection of high-resolution stellar 
spectra, our method is well-suited to complement Gaia. 


1 INTRODUCTION 


Knowing the distances of stars is fundamental for many subjects in 
astrophysics. Indeed, distances are required to understand the evo¬ 
lutionary status of stars and plac e them onto the Hertzprung-Russell 
(HR) diagram (e.g. llbenll96 7a b; Cass isil20 1 4) . which relates their 
effective temperature with their intrinsic luminosity. Although the 
effective temperature is easily determined by several observational 
techniques from spectroscopy (see e.g. iGravll 19921, for a rev iew) or 
from photometry (e.g. lAlonso. Arribas & Martmez-Rogej Il999l) . 
the intrinsic luminosity is related to the flux emitted by the star, 
which is inescapably coupled with the distance. Precisely know¬ 
ing the distance and temperature of a star allows us to separate its 
space motions, estimate physical properties such as radius, mass, 
age, and to infer internal nuclear reaction yields. Combining this 
information for many stars allows us to understand the structure 
and evolution of the Milky Way. 

The HR diagram has been succes sfully used for cluster stars 
(e.g. iGratton et al'll2003l ; |c assisil[2014h because these stars are be¬ 
lieved to have the same age, chemical composition and distance, 
and relations between their magnitudes and luminosities can be em¬ 
ployed. For individual stars in the Galactic field, the HR diagram 
can be utilised with confidence only for stars for whi ch a parallax 
has been measured dGroenewegenll2008l ; ICassis il2014l) . These stars 
correspond to the closest and brightest ones in the solar vicinity ob¬ 
served by the Hip parcos mission (jPerrv man et alJI 1997h . Very soon 
the Gaia mission l lPerrvman et al . 200lh will provide parallaxes and 
proper motions with unprecedented accuracy for one billion stars in 
the Galaxy, allowing us to construct HR diagrams beyond the solar 
vicinity. 

The fields of stellar and Galactic astrophysics are converg¬ 
ing into one field thanks to the Gaia mission and its complemen¬ 
tary groun d-based high-resoluti on spectros copic surveys, such as 
Gaia-ESO dGilmore et al.ll2012l) . 4 MOST dde Jonel 1201 ill. RA VE 
dSteinmetz et alJl2006t) . APOGEE dAllende Prieto et al.l l2008h or 
GALAH dDe Silva et alj 120151) . In this pre-Gaia era, methods 
to determine the distance of typical Galactic field stars beyond 


Hipparcos are constantly under development (jlvezic et al. 

200q; 

Burnett & Binney 2010; Binnev et al. 2014; Rodrigues et al. 

20l4 

Xue et al. 20l4l). Most methods, commonlv also referred 

to as 


“spectroscopic parallaxes", are indirect as they rely on parameters 
determined from colours or from spectra, and of course on stellar 
evolution models, to predict the position of stars in the HR dia¬ 
gram. The estimated intrinsic luminosity of a star is used to es¬ 
timate the distance, given its apparent brightness on the sky. The 
power of these methods is that they provide distances estimates for 
large numbers of field stars. Nonetheless, they require assumptions 
about extinction along the line-of-sight, but more seriously they are 
limited by systematic uncertainties in stellar evolution modelling 
and spectroscopy, which can result in distance errors of the order 
of up to 70% for metal-poor giants dSerenelli et al.l2013l) and many 
other spectral types. 

Other attempts to determine distances with better precision 
use particular types of stars with well defined intrinsic luminosities 
such as variable stars, blue horizontal branch stars or red clum p 
stars l lGroenewegen. Udalski & Bonoj l2008l ; Nidever et ahl l20l4) . 
but the luminosity of those stars is also subject to calibrations 
from stellar evolution models. An alternative idea is to compare 
the spectra of stars with those for which the parallax is known. 
ISoubiran. Bienavme & Sieberti d2003h introduced this concept to 
determine the atmospheric parameters and absolute magnitude of 
a star by computing a likelihood function of the difference between 
its spectrum and a grid of observed spectra of stars with parallaxes. 
The reported distances errors are less tha n 30% for w hen compared 
with Hipparcos parallaxes of FGK stars dkatz et al] 201 lh and less 
than 12% when red clump stars are considered I Soubiran et all 
120081) . This compares well with other methods but in this case the 
distances rely on interpolations in the parameter space to compute 
the likelihood function and not on stellar evolution and atmosphere 
models. 

We report another perspective of the concept of direct compar¬ 
ison of spectra used in ISoubiran, Bienavme & Sieberti d2003l) . but 
here we focus only on twin stars. We call this method “the twin 
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method" and distances (parallaxes) determined with it are referred 
to as “twin distances" (twin parallaxes). The method is based on 
two principles: 

• Twin stars are proposed to have the same luminosity. We iden¬ 
tify such stars by comparing their spectra, which must be identical 
according to a criterion described in Sect. [3]. 

• Because of our proposition, the difference of the apparent 
magnitudes of twin stars is directly related to the difference in 
their distances. Thus, knowing the distance of the first star from its 
parallax, allows us to directly ascend on the cosmic ladder. 

For details on the mathematical formalism of this method we 
refer the reader to the Appendix The twin method is differ¬ 
ent from other methods to determine stellar distances, in particular 
from the regular spectroscopic parallax method, because its foun¬ 
dations: First, we do not require a spectral parameter determination 
method to analyse a star spectrum with the aim to obtain its ef¬ 
fective temperature, surface gravity and metallicity. Spectroscopic 
parallax methods need this information to place the star in the FIR 
diagram which, with the help to an evolutionary model, a luminos¬ 
ity can be assigned to the star. In fact, we do not use isochrones or 
the HR diagram to determine the distances. Second, in our method 
we do not determine the value of the absolute magnitude (or dis¬ 
tance modulus) of individual stars, instead we use the relative mag¬ 
nitudes to work out the distance directly. Third, when using twin 
stars, the difference of their de-reddening is directly related to the 
observed colours of the twins, which is useful because we do not 
need to use extinction maps to correct reddened magnitudes. Fur¬ 
thermore, we consider several photometric bands in our distance 
calculation, allowing us to treat the effects of dust extinction more 
robustly as well as to assess cases with bad photometry in a given 
band. 

Although very soon Gaia will provide accurate parallaxes 
for the majority of the stars observed with high-resolution by 
current instruments, we fully agree with the point made by 
iBumett & Binnevl d2010l) . that other distance estimation methods 
complementing these parallaxes are crucial, especially for more 
distant stars. Moreover, it will remain important to have alterna¬ 
tive distance measurements of stars that are independent from par¬ 
allaxes or from model-dependent methods. This is especially true 
for stars located in remote open clusters, for which parallax mea¬ 
surements will remain challenging even for Gaia, making isochrone 
fitting to their colour-magnitude diagrams one of the only available 
methods to determine their distances. 

In this article we introduce the twin method and analyse its 
performance in a sample of FGK stars with Hipparcos parallaxes 
as well as in open clusters which have spectra taken with high reso¬ 
lution. Applications to lower resolution spectra and to stars beyond 
the reach of Hipparcos will be presented in future works on this 
subject. 


2 SPECTROSCOPIC AND PHOTOMETRIC SAMPLE 

To test the twin method we require stars to comply the following 
criteria: 

• Parallaxes from Hipparcos with accuracies of better than 10%. 

• Stellar parameters of typical FGK Milky Way stars. For that 
we took the parametersJrom the literature as compiled in the PAS¬ 
TEL catalogue dSoubiran et aljEoIoi) . For main-sequence stars we 
selected from PASTEL all stars that have temperatures in the range 
[5000, 7000] K and surface gravities in the range [3.5, 5.0] dex. For 


giants we selected from the PASTEL database stars in temperature 
range [3500,5000] K and surface gravity range [0.5, 2.0] dex. 

• They are not classified as variables by the SIMBAD database 
dWenger et ak|l2000i ). 

• The stars are not SB2 binaries, i.e. no double spectroscopic 
lines in the analyses reported by the literature compilation of PAS¬ 
TEL. 

• They have photometry in five photometric bands. 

• They have been observed with the HARPS spectrograph. 

This gives us a total of 536 field stars, whose positions in the 
sky and in the HR diagram can be seen in Fig. |T| The left panel 
shows their distribution on the sky, represented with dots. The tri¬ 
angles correspond to the open clusters initially selected for our 
study (see below). The stars are well distributed in the Southern 
Hemisphere. The right panel shows the HR diagram of our sample. 
For this illustration, the absolute mag nitude was calculated using 
the parallaxes of [van Leeuwenl d2007l) and the V magnitudes and 
B — V colours are corrected by dust extinction using the values 
provided in lESAl J1 9971 ). The gap between the main-sequence and 
the red-giant branch is a purely selection bias since we selected 
main-sequence and giants separately using temperature and logg 
cuts (see above). 

2.1 Spectra 

The HARPS instrument is fibre-fed by the Cas segrain focus of the 
3.6m telescope in La Silla lMavor et al .1 d2003l) . The spectra were 
reduced by the HARPS Data Reduction Software (version 3.1) and 
were downloaded from the advanced data product from the ESO 
archives. We chose this instrument because of the very high signal- 
to-noise ratio (SNR), high resolution and long wavelength coverage 
that a typical bright Hipparcos star spectrum has when observed 
with HARPS. 

The fact that the majority of the stars in our sample are on the 
main-sequence (see HR diagram of Fig.Q} is an observational bias, 
due to the many projects looking for planets in solar analog spectra 
in the HARPS database. Our sample includes the solar spectrum 
from the Ceres asteroid which has a signal-to-noise (SNR) of ~ 
250 for calibration purposes of our method to search for twins, i.e. 
we used it to see if we could find the solar twins in our sample that 
have been studied in other works. This spectrum was taken from 
the library of Gaia FGK benchmark star^j (see I Jofre et alj|2014al . 
for a summary), which was processed in the same way as the rest 
of the data of this study (see Sect. G.H . Our final sample contains 
536 FGK stars. 

Finally, we col lected spe ctra fro m stars belonging to 11 open 
clusters analysed bv lBlanco-Cuaresma et al.l d2015i) . The selection 
criteria for cluster stars was different than for the field stars, in the 
sense that no colour cut or another selection on the atmospheric 
parameters was done. We queried all spectra for cluster stars in 
the HARPS archives, with the only restriction to have SNR of at 
least 10. For cluster stars we stacked the spectra of the same star to 
increase the SNR. Our final sample includes 177 open cluster stars. 

2.2 Photometry 

To determine the distance from the apparent magnitudes of stars it 
is important to have as many photometric bands as possible. This is 

1 http://www.blancocuaresma.eom/s/benchmarkstars/ 
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Figure 1. Our sample of Hipparcos stars with HARPS spectra. A: Distribution of stars across the observable Southern sky from La Silla Observatory, Chile. 
Triangles correspond to positions of clusters. The sinusoidal shape drawn by the open clusters reflects the location of the Galactic disk. B: HR diagram of the 
Hipparcos sample. The absolute magnitudes were computed using the provided parallax and the apparent magnitudes B and V are in the Johnson-Cousins 
filter system. The main-sequence and red giant branch are visible. 


crucial in our analysis, as the distance is obtained from the differ¬ 
ence of magnitudes, which needs to yield consistent results when 
different photometric bands are taken into account. Thus, when the 
results obtained with different photometric bands agree, we ensure 
the reliability of the distance obtained, as well as the photometric 
quality of the star. Another very important aspect of having sev¬ 
eral photometric bands for our analysis is the effect of extinction 
or saturation, which is different for the different bands. When de¬ 
termining distances from different photometric bands, the effects 
of extinction or saturation are less dependent on our final results. 
Furthermore, it is possible that some stars are missclassified vari¬ 
able stars. Different photometric measurements are useful for di¬ 
agnosing if this is the case. For these reasons, the determination 
of dista nce was done using the photometric ban ds from H ippar- 
cos Hp ( Ivan Leeuwenl20070 . Johnson-Cousin B, V dESAl 1997D and 
2MASS J, Ks dSkrutskie et al .1120061) . 


The ratio of total-to-selective extinction R(A') was assumed to 
have five diffe rent values corresp ondin g to the five different band s 
as taken from iFitzpatrickl d 19991) and Yuan. Liu & Xiand d2013l) . 
The values are R (H p ) = 3.1, R(-B) = 3.8, R(U) = 3.1, 
R(J) = 0.7 and R(/Ts) = 0.3. Note that the reddest Ks band 
has a very small extinction, with R (Ks) of 0.3, while the bluest 
B band has a significantly larger extinction, with R (B) of 3.8. It 
is important to remark that our present analysis for field stars is 
restricted to Hipparcos stars, which are close by and are not sig¬ 
nificantly affected by extinction dLallement et al.ll2003l) . The stars 
of our sample are mostly nearby disk stars, and the assumption of 
constant R seem to be reasonable in the five photometric bands 
considered. We will see in Sect. l4.2l that using five bands also help 
us to treat extinction properly in open clusters. Hence, we prefer to 
use our generic form consistently throughout our work by applying 
the extinction correction factor of Eq. lAbl in our analysis. 


3 PROCEDURE TO FIND SPECTROSCOPIC TWINS 
3.1 Spectral pre-processing 

The spectra needed to be pre-processed before the analysis. We 
carefully normalised the spectra, corrected for the radial veloci¬ 
ties of each star, removed telluric absorption, and resampled the 
spectr a to a common wav elength range using the tools described 
in lBlanco-Cuaresma et all J20l4hl) . This was done homogeneously 
in order to ensure that the relative comparison of the spectra is as 
objective as possible. 

The normalisation was done by fitting splines of second or¬ 
der locally to the spectra as described in IBlanco-Cuaresma et alj 
d2014bl) . The radial velocity correction was done by cross- 
correlating the spectra against a template rest-frame. This template 
was either an atlas of the spectrum of the Sun or of Arcturus, de¬ 
pending on the colour of the star. If the star had B — V < 1.0 mag, 
the star is most probably a dwarf and we used the template of the 
Sun. Otherwise if the star had B — V > 1.0 mag, then the star is 
probably a giant and the template of Arcturus was more suitable 
for the cross-correlation. Both atl ases a re included with the tools to 
prepare the library of lBlanco-Cuaresma et al J d20 14bl) . 


3.2 Spectral analysis 

Following previous spectroscopi c studi es of solar twins (e.g . 
I Melendez, Dodds-Eden & Roblesl 120061 : iRamfrez et al] I20l4 
Datson, Flynn & Portinari 2014) we measured the strength of 
atomic absorption lines, i.e. we worked with equivalent widths 
(EWs) by differentiating each and studied them as a function 
of excitation potential (E.P.). To account for broadening effects 
caused by e.g. rotation, we also investigated the depth of each 
line by calculating the value of their minima with respect to 
fluxes at ±0 .02 nmjrom the centre of the line, similar to what 
IMelendez. Dodds-Eden & Roblesl te.g. l2006 ) have done. 

We employed the atomic data for Fe, Si, Mg, Ca, Ti, Sc, V, 
Cr, Mn, Co a nd Ni used for the Gaia benchmark stars djofre et alJ 
I2014bl . l2015l) . which provide a substantial number of transitions in 
FGK stars. We calculated EWs of 423 atomic lines listed in these 
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Figure 2. Example of equivalent width differences of four pair of stars as a function of excitation potential. A linear regression fit is performed to the data of 
neutral lines (blue open circles) and ionised lines (red filled circles), separately. The slope and zero point of the fit is indicated in each case, with the error of 
the slope in parenthesis. The names of the stars from the pair are indicated in the right hand side of each panel. Two examples are shown, where panels Al and 
A2 show regular stars and panels B1 and B2 show twin stars. 


papers, which have heen carefully selected to be unblended lines 
for typical FGK Milky Way stars. To assess the systematic error of 
the EW measurement, we employed two different codes to measure 
the EWs: iSpec iBlanco-Cuaresma et alj2014ah and SMH Jcasevl 
l20l4 ). We selected the lines with EWs between 15 and 120 mA 
showing a good agreement between both methods. These lines are 
strong enough to have accurate EW measurements in both methods, 
but weak enough for the line not to saturate. The EWs were then 
compared between two spectra, which reduced our data size from 
262,000 pixels to <423 points (per star), facilitating a faster and 
focussed comparison. 

We analysed the trend and offset of the differences of EWs as 
a function of E.P. for each combination of spectra, which was done 
by fitting a line to the data by linear regression. To obtain a robust 


estimate of slope of our linear fit, we performed bootstrapping by 
fitting a dataset of randomly selected 80% of lines 1000 times. The 
final slope was the mean of the 1000 fits obtained, with its standard 
deviation as its uncertainty. 

Twin stars are required to have fits with zero slope and zero 
offset (value of the fit at E.P = 0) of both, ionised and neutral lines. 
Since the linear fits have a dispersion around zero, we carefully 
studied the values of slopes and offsets for our complete dataset. 
After trying several different combinations and performing visual 
inspection of the selected twin pair spectra, we found that a good 
compromise of spectroscopic twins and number of selected pairs 
was obtained under the following criteria: 
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Wavelength [nm] 

Figure 3. Example of twins observed with HARPS for stars with FGK spec¬ 
tral classes. Coloured continuous lines represent one star, and black dotted 
lines represent the twin. Effective temperature decreases from top to bot¬ 
tom. 

(i) The slope and offset do not differ by more than 0.3 mA/eV 
and 2.5 mA for neutral lines, respectively. 

(ii) Error of the slope of neutral lines of less than 0.1 mA/eV. 

(iii) The offset between the ionised and neutral lines is not more 
than 2.5 mA for ionised lines with respect to neutral lines. 

(iv) Mean difference of line depths of less than 10%. 

Examples of how EWs compare for different stars are found 
in Fig. [2] The open blue circles represent neutral lines and the filled 
red circles represent ionised lines. The regression linear fits to the 
data are represented with blue and red dashed lines, for the neutral 
and ionised lines, respectively. The values of the slope and offset 
are indicated, with the error of the slope in parentheses. The first 
two panels show the star HIP 36512 compared to a regular star 
(Panel Al) and a twin star (Panel Bl), respectively. The last two 
panels show another example, but for the star HIP 41479, in which 
Panel A2 and B2 show a regular and a twin star, respectively. For 
the regular cases (Panels Al and A2), we see clear trends and off¬ 
sets in the linear fits, while for the twin stars (Panels Bl and B2), 
the fits follow the zero line, as defined with our criteria above. 

Three examples of twin pairs identified using our method can 
be seen in Fig. [3] Each example shows a spectral region around the 
Mg I b triplet, a feature known to be very sensitive to stellar param¬ 
eters. For this illustration, one star of each twin pair is plotted with 
a continuous coloured line, and the other one is plotted with black 
dotted line. The spectra are indistinguishable within the errors; 
without knowing the coordinates and magnitudes, it is impossible 
to disentangle if both spectra correspond to two independent stars 
or the same one observed in different conditions. 

A total of 175 twin pairs were found using these criteria. 
These pairs include triplets or quadruplets, partly because of the 
biases in the HARPS sample towards solar analogs. Although this 
number seems very large, it is small given the total amount of pairs 
analysed in this work (nx(n+l)/2 where n is the number of stars 
in the dataset.) In this particular sample, we find 175 twin pairs out 
of a total of 12844 comparisons, i.e. 1.36% of the total pairs re¬ 


HR diagram 



Figure 4. HR diagram of the Hipparcos sample. Open circles: the entire 
sample of 536 stars, filled circles: stars with twins. Blue colour corresponds 
to spectroscopic twins while red colour indicates the stars for which a dis¬ 
tance could be determined with a good consistency over several photometric 
bands. 

suited in twin stars. Their location in the HR diagram can be seen 
in Fig. [4] The filled circles (blue and red) are stars with twins in 
our sample that satisfy the spectroscopic criteria. The red circles 
correspond to our final set of twins, which had good photometry 
(see Sect. ns. From Fig.[4]we can see that most of the twins found 
are well distributed in the main-sequence, although some of them 
are subgiants and red clump stars. We could not find twins for the 
cool giants of our sample, which is not surprising when seeing how 
they distribute in the HR diagram. One can find few stars with sim¬ 
ilar colour and luminosity at these temperatures. The few spectro¬ 
scopic twins found on the red giant branch had photometry in dif¬ 
ferent bands that was not consistent enough for reliable distance 
measurements. 


4 RESULTS ON DISTANCES 

4.1 Twin parallaxes within the Hipparcos sample 

We considered five different photometric bands to determine the 
distance from Eq. ( lA6l for each band. Only the values having 
a standard deviation smaller than 20% of the measured parallax 
due to different photometric bands were considered to have accu¬ 
rate photometry, thereby assuring reliable results. These stars are 
marked with red colour in Fig. [4] One can see that very few spec¬ 
troscopic twins yielded inconsistent distances when considering the 
five different photometric bands. Uncertainties of the measured par¬ 
allax of the reference star, as well as the errors in the photometry 
were considered for calculating the final uncertainty in our paral¬ 
laxes. This was done using Monte Carlo simulations on the errors 
in the measured parallax of the reference star and the photometry. 
We found that the photometric error contributed to negligible un¬ 
certainties with respect to the initial parallax error and the standard 
deviation of the five different photometric bands. 

The mean uncertainty in the Hipparcos parallaxes of our 
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COSMIC LADDER OF HIPPARCOS TWINS 
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Figure 5. Difference of apparent magnitude of twin stars respect to the differences in their Hipparcos parallaxes. Each panel shows the difference for different 
photometric bands. The red color correspond to the stars when corrected by extinction. The blue solid line corresponds to the line with slope of 5. The figure 
shows that Eq, IA6l is satisfied for twin stars, i.e, the proposition they have the same luminosity is correct. 


sample is 2.75%, which propagated to a mean uncertainty in 
the twin parallaxes of 3.21%. The mean error due the different 
photometric bands is of 1.72%, which gives us a mean total 
uncertainty of 3.46% in the twin parallaxes. 

For the Hipparcos sample, 175 twin pairs with reliable pho¬ 
tometry were found. We show in Fig. [5] that our proposition that 
spectroscopic twin stars have the same luminosity and Eq. lA6l can 
be used to determine distances. The figure shows the logarithm of 
the ratio of the twins’ parallaxes as measured by Hipparcos, versus 
the difference in the apparent magnitudes. The solid line of slope 5 
is overplotted with blue. Each figure uses the difference of magni¬ 
tudes in the five different bands considered in this work. The data 
has been plotted twice in each panel, the first time we do not con¬ 
sider extinction (essentially Eg I A1 b and the data is plotted in black 
crosses, while the second one we consider extinction (essentially 
Eg I A6b with red filled circles. 

Our selected spectroscopic twins satisfy our relation very well, 
showing like that the difference in the apparent magnitudes is di¬ 
rectly related to the difference in their distances, ascending like this 
the cosmic ladder. The figure also tell us that the the corrections 
due to extinction are normally very small. The corrected magni¬ 
tudes (red filled circles) have a scatter along the line of slope 5 that 
agrees over a range of 3 orders of magnitude with respect to the 
uncorrected magnitudes. Also, the scatter is very similar within 3 
orders of magnitude for all bands. This is expected as we selected 
our twins based on stars giving consistent distances between all 
bands. We comment here that although we show that extinction is 
in this case negligible, when applying the method to datasets that 
expand beyond the local bubble this will not be the case, the general 
form of Eq. lA6l should be used. __ 

The comparison of the Hipparcos parallax l lvan Leeuwenl 
120071 ) with our results (twin parallax) are displayed in Fig. [6] In 
Panel A we show the direct comparison of the Hipparcos parallax 
and the twin parallax for the sample. The distribution of the de- 
viatior@ between our results and Hipparcos is shown in Panel B. 

2 We define deviation as dev(plx) = (1 — ntwin/^hip) X 100 with 
ntwin. the twin parallax and n^p the Hipparcos measurement. 


Assuming this distribution can be represented by a Gaussian, its 
mean and standard deviation are indicated at the top of Panel A. 
The agreement is excellent with a standard deviation of of 7.5%. 

In Panel C of Fig. [6] we display the deviation between Hip¬ 
parcos parallaxes and our results as a function of distance, while 
in Panel B we display the deviation as a function of colour. In 
both panels the dashed lines correspond to the standard deviation of 
7.5%, while the dotted lines represent the deviation of 20%. We can 
see that all the determined twins parallaxes in our sample deviate 
from the Hipparcos ones by less than 25%. We can also see that, 
although our sample is focused on stars that are very close by, the 
agreement between our results and those of Hipparcos is indepen¬ 
dent on the distance. Similarly, we determine accurate distances for 
both dwarf and giant stars, i.e . for all colours in our sample. 

Recently. [Santiago et al.ld2015 ) presented a spectroscopic par¬ 
allax method to determine distance in which they analyse a sample 
of stars that overlap with our sample. They took approximately 500 
main-sequence stars observed with HARPS considering the stel¬ 
lar parameters derived by ISousa et~ahl <1201 ll) and obtained agree¬ 
ments with Hipparcos parallaxes of the order of 20%. Our results, 
although providing distances for a subsample, agrees better with 
7.5%. 


4.2 Application to open clusters 

The stars in open clusters do not have parallax measurements, 
so our results would mostly compare with the model-dependent 
method of isochrone fitting to the main-sequence of the open clus¬ 
ter. 

To search for twins in open clusters we had to neglect the sec¬ 
ond criterion listed in 13.21 The reason is that many of the spectra 
in the clusters had significantly lower SNR (SNR ~ 15) than the 
Hipparcos stars, so the measurement of EWs was more uncertain, 
which produced a larger scatter in fitting the line to the differences 
in EWs as a function of E.P. We visually inspected the spectra 
of each of the twin candidates found using the other three crite¬ 
ria listed in !3.2l The initial sample of twins found in open clusters 
was reduced with an extra requirement of stars having 2MASS and 
Johnson-Cousins photometry. We clarify this point here as not all 
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Figure 6. Distances obtained with the twin method compared to Hipparcos measurements. A: The one-to-one relation is shown with a continuous line and the 
parallaxes of the twins are shown as dots. B: The distribution of deviations of the parallaxes. C: Deviation of the distance is plotted as a function of distance 
with dashed line representing 1 cr of the deviation distribution of panel B. Dotted line represents the limit of 20% deviation. D: as panel C with distances plotted 
as a function of B — V colour. 


cluster stars initially selected from iBlanco-Cuaresma et al . < 120151 ) 
had photometry in the four bands required for our present study, in 
which we need to compare several photometric bands in the same 
way as the Hipparcos sample. This requirement is especially impor¬ 
tant in clusters, where extinction and binary fraction can become 
significantly larger than in the field. Only stars for which we ob¬ 
tained a better agreement than 20% between the distances obtained 
using the four different photometric bands were selected (cluster 
stars do not have H p photometry). The clusters M67 and NGC2360 
had sufficient member twins of our Hipparcos sample with the re¬ 
quired photometric accuracy. Thus, only they were suitable for us to 
provide a robust value for their distance. For the rest of the clusters, 
we either did not find any twin with our library, or we found very 
few, but most of them lacking 2MASS photometry or not having 
the radial velocities of the other cluster members. 

It is important to comment that most of the cluster spectra cor¬ 
respond to giants in the cluster, and our Hipparcos database con¬ 
tains mainly dwarfs. Enlarging the overlap between field and clus¬ 
ter stars can not be easily done since this study was limited to what 
is available in the public HARPS dataset. On one hand, most of 
the field stars are dwarfs because of the intensive search of exo¬ 
planets around solar-like stars. On the other hand, because clusters 
are more distant, giant stars are bright enough for HARPS to get 
high SNR. Finding more than 5 twins with good photometry in two 
clusters using only archive data is indeed very encouraging. The 
member stars with twins of Hipparcos are plotted with red colour 
in the colour-m agnitude diagra ms of Fig.Qas taken from WFBD/'Q 
combined with lPasquini et all d2008l) for the twins found in M67. 

From the diagram we can see that the twins found in 
NGC2360 are red clump giants while the twins found in M67 are 


3 http://www.univie.ac.at/webda/ 


NGC2360 M67 




Figure 7. Colour-magnitude diagram of the two clusters in this study. The 
filled red circles correspond to the star clusters that have a twin in the Hip¬ 
parcos catalogue. For NGC2360 we only have twins for giant stars, while 
for M67 we found giants and dwarfs. 

both, main sequence stars and giants. There is in fact only one gi¬ 
ant in M67 with an Hipparcos twin, which is not at the red clump. 
The distances obtained for each of the stars are displayed in Fig. [8] 
as a function of colour. The error bars are the square root of the 
quadratic sum of the error due to the uncertainty of the parallax 
from the Hipparcos twin and the error due to the distance obtained 
from the different photometric bands. 

The solid line in Fig. [8] represents the median of 
the distances obtained for each star individually, the 
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Figure 8. Distance of the cluster M67 and NGC2360 in Panel A and B, 
respectively. The dashed horizontal lines represent values reported in the 
literature and continuous lines represents our final distance result. Dots cor¬ 
respond to the distances obtained for each twin found in the cluster, as a 
function of colour. The error bars represent the standard deviation of the 
results obtained from different photometric bands. 


dashed lines represent the values reported in the literature 
l ^varoe^shmai^fcAnthon^Xvaro^. j_997|^_ Kharchenko & et al\ 
120131 : IGiines. Karatas & Bonattol l2012t IPasauini et al.l 120081) . 
Not all literature studi es liste d here report distances for both 
clusters. IPasauini tttak _i 2008 )_ determi ned the distance only 
for M67 while iGiines Karatas & Bonattol l2012h only for 


NGC2360. In Fig. [8] we show in the legend the agreement with 
Twarog. Ashman & Anthony-Twarogj (l997j), who analysed both 
clusters and agrees with us better than 3%. Note that all literature 
values except that one of IPasauini et al.1 J 2 OO 8 I) determined the 
distance using evolutionary models by fitting the main-sequence to 
the colour-magnitude diagram of the cluster. Our method is more 
similar to IPasauini et al.l d2008l ). which considered the averaged 
magnitude of ten solar analogs in M67 and the absolute magnitude 
of the Sun to compute the distance modulus. The agreement with 
Pasquini et al] <20081) is remarkable. Also note that the discrepancy 
in distances for NG C2360 using main-se q uence fitting exceeds 
35% for the case of I Kharchenko & et a 3 d20l3l) respect to the 
other works. 

It is curious that our result obtained from the reddest (giant) 
star in M67 is different than those obtained from the bluer 
(main-sequence) stars. Although we did not find an offset for giant 
stars in our Hipparcos sample, measuring distances to red giants 
can b e very challenging (e.g. ISerenelli et al.l 120131 : iBinnev et all 
120141) . The shape of the red giant branch (Fig. [7J, demonstrates 
how a slight change in colour (temperature) implies a large change 
in magnitude (luminosity). Having a larger sample of giants is 
needed to understand why this star gave a lower distance than its 
dwarfs siblings in M67. Such insights could be revealed from a 


stellar survey like APOGEE or Gaia-ESO, where large samples of 
giants are observed in clusters and in the field. 


5 DISCUSSION AND CONCLUSIONS 

The twin method is simple, and is based on two fundamental as- 
trophysical principles: the apparent magnitude of a star decreases 
with distance; and if two stars are equivalent then their spectra are 
identical. When several photometric bands give consistent magni¬ 
tudes and the parallax of the reference star is well measured, the 
twin method has important applications in stellar and Galactic as¬ 
trophysics beyond the direct determination of distances. 

It becomes a powerful method to complement future Gaia par¬ 
allaxes. It also allows in-depth studies into astrophysical mecha¬ 
nisms affecting stellar luminosities. For example, Fig.[6]shows that 
there are cases with distances that deviate more than 15% from 
Hipparcos. Our twins have been found under the assumption that 
the flux around the spectral lines of 11 chemical elements agree. It 
is possible that the absorption produced by other elements, such as 
carbon or neutron-capture elements, might change opacities such 
that the evolution, and thus the luminosity, is slightly affected. In 
our procedure to find twins, every EWs has the same weight when 
we perform the fitting of the trends. Other elements that are not in¬ 
cluded in the selecti on of lines of the Gaia FGK benchmark stars 
in jjofre et alJl2014bL 2015) are less numerous in the optical spec¬ 
tra. Including the EWs of some lines from these elements would 
produce an effectively unaffected fit of the trends. The_study of the 
effect of variation of other elements not listed in Jjofre et al.l2014bl . 
2015) in the deviation of the distance goes beyond the purpose of 
the current paper, which consists in introducing the concept of us¬ 
ing twin stars for the determination of distances. We have shown 
here that by comparing equivalent widths of the elements Mg, Si, 
Ca, Ti, Sc, V, Cr, Mn, Fe, Co and Ni we are able to find stars that are 
equal enough to determine distances within 7.5% agreement with 
Hipparcos parallaxes. 

The twin method also provides an attractive approach to quan¬ 
tify the effects of dust extinction in the Milky Way. The fac¬ 
tor R(X) used in Eq. j A6b is assumed to be constant in the 
disk, but this value origin ates from calibrations lFitzpatrickl ll999l : 
lYuan. Liu & Xianel 120131) . One can refine these calibrations by 
identifying twins with good photometry and parallaxes that are lo¬ 
cated in different parts of the Galaxy. Furthermore, one can test if 
R(.Y) is still constant outside the disk. Such analysis could be done 
when Gaia releases its parallaxes, giving us accurate distances for 
stars beyond the solar vicinity and the disk. 

It is worth to comment that our method, by being restricted to 
twin stars, provides distances to a much smaller number of stars. 
This happens however to several other methods that are also re¬ 
stricted to particular kind of stars, such as distances to variable 
stars (e.g^Genovali_etal 2013j) or distances using asteroseismic 
data dSilva Aguirre et aHl2013l) . In this era of large datasets of stel¬ 
lar spectra and with Gaia parallaxes coming soon, samples of twin 
pairs with one star having a parallax will become numerous in the 
near future. 

With the relatively small sample of 536 high resolution spec¬ 
tra of FGK stars we found 175 twin pairs. We showed that with 
accurate photometry and parallaxes we can retrieve twin parallaxes 
with typical errors of 3.25%, which is of the order of the Hippar¬ 
cos uncertainties for our sample. Our results agree within 7.5% for 
most of the cases with the measurements of Hipparcos. From this 
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sample we also determined the distances for the clusters M67 and 
NGC2360, finding agreement of ~ 3% w it h sta ndard techniques 
such as that one of tTwarog. Ashman & Anthonv-Twarogi Il99 7i). 
Our main result can be seen in Fig. [5] which demonstrates that our 
basic proposition, that twin stars have the same luminosity (i.e. ab¬ 
solute magnitude) is verified. 

Currently, high-resolution spectra of Hipparcos stars can be 
found in several public archives. Extending our sample of refer¬ 
ence stars will allow us to find more twins of Hipparcos in on¬ 
going high-resolution spectroscopic surveys such as Gaia-ESO, 
APOGEE or GALAH, and determine distances of fainter objects 
that do not benefit from direct parallaxes. With Gaia measuring 
the most accurate parallaxes in history, millions of stars already 
observed with high-resolution spectra will soon have accurate dis¬ 
tances, allowing us to climb the cosmic distance ladder. Indeed, 
with the next generation of large telescopes, such as E-ELT, high 
resolution spectra of objects beyond the reach of Gaia by several 
magnitudes, will become available. We look forward to applying 
the twin method to stars too distant for accurate parallaxes from 
Gaia, but for which good spectra can be obtained with large tele¬ 
scopes, and helping to understand fundamental astrophysical pro¬ 
cesses from the deepest parts of our Galaxy. 

Acknowledgments 

It is our pleasure to thank T. Masseron, K. Hawkins and J. Pen- 
ner for fruitful discussions on the subject. The authors thank 
Chris Flynn for the good suggestions made in the referee report, 
which helped to improve the paper significantly. This work was 
partly supported by the European Union FP7 programme through 
ERC grant number 320360. C. W. acknowledges Leverhulme Trust 
through grant RPG-2012-541. Based on data obtained from the 
ESO Science Archive Facility. This research has made use of the 
SIMBAD and WEB DA database. P. J. and T. M. thank Damian and 
Kian for motivating this work. 


REFERENCES 

Allende Prieto C. et al., 2008, Astronomische Nachrichten, 329, 
1018 

Alonso A., Arribas S., Martlnez-Roger C., 1999, A&AS, 140, 261 
Binney J. et al., 2014, MNRAS, 437, 351 
Blanco-Cuaresma S. et al., 2015, A&A, 577. A47 
Blanco-Cuaresma S., Soubiran C., Heiter U., Jofre P, 2014a, 
A&A, 569, Al 11 

Blanco-Cuaresma S., Soubiran C., Jofre P, Heiter U., 2014b, 
A&A, 566, A98 

Burnett B„ Binney J., 2010, MNRAS, 407, 339 
Casey A. R.. 2014, ArXiv e-prints 

Cassisi S., 2014, in EAS Publications Series, Vol. 65, , pp. 17-74 
Datson J., Flynn C., Portinari L„ 2014, MNRAS, 439, 1028 
de Jong R., 2011, The Messenger, 145, 14 
De Pascale M., Worley C. C., de Laverny P, Recio-Blanco A., 
Hill V., Bijaoui A„ 2014, A&A, 570, A68 
De Silva G. M., Freeman K. C., Bland-Hawthorn J., Martell S., de 
Boer E. W.. Asplund M., Keller S., 2015, ArXiv e-prints 
ESA, ed„ 1997, ESA Special Publication, Vol. 1200, The HIP¬ 
PARCOS and TYCHO catalogues. Astrometric and photometric 
star catalogues derived from the ESA HIPPARCOS Space As¬ 
trometry Mission 

FitzpatrickE. L„ 1999, PASP, 111, 63 


Genovali K. et al., 2013, A&A, 554, A132 
Gilmore G., Randich S., Asplund M., Binney J., Bonifacio P., 
Drew J., Feltzing S., Ferguson A., 2012, The Messenger, 147, 
25 

Gratton R. G., Bragaglia A., Carretta E., Clementini G., Desidera 
S., Grundahl F„ Lucatello S„ 2003, A&A, 408, 529 
Gray D. F., 1992, The observation and analysis of stellar photo¬ 
spheres. 

Groenewegen M. A. T., 2008, A&A, 488, 935 
Groenewegen M. A. T., Udalski A., Bono G., 2008, A&A, 481, 
441 

Giineg O., Karata§ Y., Bonatto C., 2012, New, 17, 720 
Iben, Jr. I., 1967a, Science, 155, 785 
Iben, Jr. I., 1967b, ApJ, 147, 624 
Ivezic Z. et al., 2008, ApJ, 684, 287 

Jofre P, Heiter U., Blanco-Cuaresma S., Soubiran C., 2014a, in 
Astronomical Society of India Conference Series, Vol. 11, As¬ 
tronomical Society of India Conference Series, pp. 159-166 
Jofre P. et al., 2015, ArXiv e-prints 
Jofre P. et al., 2014b, A&A, 564, A133 

Katz D., Soubiran C., Cayrel R., Barbuy B., Friel E., Bienayme 
O., Perrin M.-N., 2011, A&A, 525, A90 
Kharchenko N. V., et al., 2013, A&A, 558, A53 
Lallement R„ Welsh B. Y„ Vergely J. L„ Crifo F„ Sfeir D„ 2003, 
A&A, 411.447 

Mayor M.. Pepe F., Queloz D., Bouchy F., Rupprecht G., Lo Curto 
G., Avila G., Benz W., 2003, The Messenger, 114, 20 
Melendez J., Dodds-Eden K., Robles J. A., 2006, ApJL, 641, 
L133 

Nidever D. L., Bovy J., Bird J. C., Andrews B. H., Hayden M., 
Holtzman J., Majewski S. R., Smith V., 2014, ApJ, 796, 38 
Pasquini L., Biazzo K., Bonifacio P, Randich S., Bedin L. R., 
2008, A&A, 489, 677 

Perryman M. A. C. et al., 2001, A&A, 369, 339 

Perryman M. A. C. et al., 1997, A&A, 323, L49 

Ramirez I. et al., 2014, A&A, 572, A48 

Rodrigues T. S. et al., 2014, MNRAS, 445, 2758 

Santiago B. X. et al., 2015, ArXiv e-prints 

Schlafly E. F„ Finkbeiner D. P, 2011, ApJ, 737. 103 

Schlegel D. J., Finkbeiner D. P, Davis M., 1998, ApJ, 500, 525 

Schultheis M. et al., 2014, AJ, 148, 24 

Serenelli A. M., Bergemann M., Ruchti G., Casagrande L., 2013, 
MNRAS, 429, 3645 

Silva Aguirre V., Casagrande L., Basu S., Campante T. L., Chaplin 
W. J., Huber D., Miglio A., Serenelli A. M., 2013, Astronomis¬ 
che Nachrichten, 334, 22 

Skrutskie M. F., Cutri R. M., Stiening R., Weinberg M. D., Schnei¬ 
der S., Carpenter J. M„ 2006, AJ, 131, 1163 
Soubiran C., Bienayme O., Mishenina T. V., Kovtyukh V. V., 
2008, A&A, 480, 91 

Soubiran C., Bienayme O., Siebert A., 2003, A&A, 398, 141 
Soubiran C., Le Campion J.-F., Cayrel de Strobel G., Caillo A., 
2010, A&A, 515, All 1 

Sousa S. G., Santos N. C., Israelian G., Mayor M., Udry S., 2011, 
A&A, 533, A141 

Steinmetz M., Zwitter T., Siebert A., Watson F. G., Freeman K. C., 
Munari U„ Campbell R„ Williams M„ 2006, AJ, 132, 1645 
Twarog B. A., Ashman K. M., Anthony-Twarog B. J., 1997, AJ, 
114, 2556 

van Leeuwen F., 2007, A&A, 474, 653 
Wenger M. et al., 2000, A&AS, 143, 9 






10 P. Jofre el. al. 


Xue X.-X., Ma Z., Rix H.-W., Morrison H. L., Harding P., Beers 
T. C., Ivans 1.1., 2014, ApJ, 784, 170 
Yuan H. B„ Liu X. W„ Xiang M. S., 2013, MNRAS, 430, 2188 


APPENDIX A: FORMALISM OF DISTANCE 
DETERMINATION OF TWIN STARS 

By restricting ourself to use twin stars, we can apply the classi¬ 
cal expression described in astronomy textbooks to measure stel¬ 
lar distances directly. Suppose there are two stars in the Galaxy at 
different locations. Assuming star 1 and star 2 are twins, they have 
the same luminosity. Thus their respective apparent magnitudes mi 
and m 2 can be related to their distances d\ and d 2 via the relation 

mi — m 2 = — 5 log(di/d2) • (Al) 

Depending on the location and distance of each star in the Galaxy, 
their apparent brightnesses will suffer from interstellar extinction. 
Suppose both stars are observed with an arbitrary filter X, their 
apparent magnitudes m; ,(* — 1,2), may assumed to be 

m; = Xi — A(X)i , (A2) 

where A (X)i represents the correction due to extinction in the 
interstellar medium, may be expressed as function of the colour 
A' — Y. 


A(X) = R(X)E(X - Y) , (A3) 

where R(X) is the ratio of the of total-to-selective extinction in 
filt er A', an d E(X — Y) is the de-reddening of the stars’ colour (e.g. 
iFitzpatrickll 19991 ). The value of the ratio of the of total-to-selective 
extinction R(X) may be assumed constant in the Galactic disk (e.g. 
Fitz patrickll 19991 : 1 Yuan. Liu & Xian3l2013l : ISchlaflv & Finkbeineil 
201 If ), therefore Eq. ( I A1 b becomes 


Ai-A 2 -R(A) 


E(X-Y)i-E(X-Y) 2 


—5 log(di/d 2 ) .(A4) 


It is seen that the difference of the logarithm of distances becomes 
a function of the difference of the apparent brightnesses and the 
difference of the de-reddening of both stars. Since the de-reddening 
is a correction to the intrinsic colour (A — Y)o of the stars, and 
since the stars are twins, their intrinsic colour is the same, which 
implies 

E(X - Y)i - E(X - Y) a = (A - Y) i - (A - Y ) 2 . (A5) 

Hence, for twins stars the difference in their observed colours is 
the same as the difference in their de-reddening. This is important 
because by considering the colours we do not need to consider de¬ 
reddening factors from extinction maps or calibrations of extinction 
in different photometric bands, which would bring external uncer¬ 
tainties to our analysis. Using Eq. jA5b . Eq. dA4b becomes our final 
expression 


Ai-Aa-R(A) 


(Y-X)i 


(Y-X ) 2 


-51og(di/d 2 ) .(A6) 


This formula can be used to determine the distance of one star 
of the twin pair knowing the photometry of both of them, the value 
of R(X), and the distance of its twin counterpart. Note that the 
only dependency on extinction in this formula is the factor R(X). 
The demonstration of this formula can be seen in Fig. [5] with a 
sample of Hipparcos twin stars, in which we verify that in a log-log 
plot, the data lie on a line with slope 5. In the figure we also plot 


the relation for Ea. IAll and we see that the extinction correction is 
negligible. This is expected since the Hipparcos twin stars lie in the 
local bubble. This is not the case for stars beyond the local bubble, 
for which the correction term of Eq. lA6l should be used. 

It is worth commenting here that the factor R(X) is indi¬ 
rectly used by any other method determining distances based on 
apparent magnitudes of stars. Usually, to know the absolute mag¬ 
nitude of an individual star (or its distance modulus) would require 
the correction of the observed magnitude from extinction, which is 
done applying Eq. to the magnitudes. The factor E(A' — Y) 
is commonly taken from ex ti nction map s, such as those from 
ISchleeel. Finkbeiner & Davisl dl998h and iSchlaflv & Finkbeinei] 
( 1201 IB for the ( B — V) colour. Extinc tion maps are different from 
each other (see lSchultheis et ail 1 20141 for such analysis), as well 
as the empirical corrections applied to E(B — V ) for correcting 
extinction of magnitudes taken with other photometric filters (e.g. 

I Yuan. Liu & Xianal2013l) . Extinction coefficients c an also be de¬ 
termi ned together with the distance indirectly (e.g. Sant iago et alj 
jml). By using the differential magnitudes in different photomet¬ 
ric bands as well as the fact that twins have intrinsically the same 
absolute magnitudes and colours, our expression to determine dis¬ 
tances fEo, IA6b is independent on extinction maps or coefficients. 


APPENDIX B: STELLAR ATMOSPHERIC PARAMETERS 

As a consistency check, we compared the stellar parameters (effec¬ 
tive temperature, surface gravity, metallicity and abundance of a 
elements with respect to iron) of our twin s as obtained within the 
AMBRE project for the HARPS dataset JPe Pascale et alj|2014l ), 
for which 135 pairs had both stars with well determined parame¬ 
ters. The remaining stars either the AMBRE method to determine 
the parameters did not converge to a physical solution, or our stars 
were observed after the AMBRE project collected the data. The 
atmospheric parameters of the twin sample are shown in Fig. HU 
where in the x-axis we plot the parameters of the first star and in 
the y-axis the parameters of the twin. The agreement is excellent, 
with negligible offsets and standard deviation that are smaller than 
typical errors on the stellar parameters. From this test we note the 
following: 

• For twin stars the AMBRE parameters agree remarkably well, 
confirming the good performance of this method to determine pa¬ 
rameters. 

• Our method yields consistent results with standard methods to 
determine stellar parameters. This validates the concept that a pure 
comparison of observed spectra serves to de termi ne parameters, 
which has been already shown by ISoubiran. Bienavme & Siebert) 
d2003l) . This also validates that the 11 elements and 423 atomic 
lines chosen for the comparison are a good proxy for the determi¬ 
nation of stellar parameters. 

• Retrieving surface gravities more accurately than 0.1 dex from 
pure spectroscopy is extremely challenging. Even in this case, 
where we have high signal-to-noise and high-resolution spectra that 
are practically identical, we find differences in the AMBRE gravi¬ 
ties of ~ 0.3 dex or more. 

• According to AMBRE parameters, our twin stars span a wide 
range of temperatures and metallicities, although they are all con¬ 
centrated towards high gravities. The reason is probably that our 
initial sample contains mostly main-sequence stars (see Fig. Q}- 
Nonetheless, from the HR diagram in Fig. [4] we see that few gi¬ 
ant twins were found, but unfortunately AMBRE parameters for 
them were unavailable. 
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Figure Bl. Comparison of effective temperatures, surface gravities, metallicities and a abundances for the twin pairs found in the HARPS sample. The 
parameters were determined homogeneously and automatically as part of the AMBRE project. 


We state here that our analysis based on EWs is not intended to de¬ 
rive atmospheric parameters. The quantities we look at, such as the 
slope of neutral lines EWs or the difference between ionised and 
neutral lines EWs, are related to the differences in the atmospheric 
parameters. We know that atmospheric parameters are crucial for 
distinguishing the different kind of stars from the spectra. More¬ 
over, from e.g. regular spectroscopic parallax methods, we also 
know that these parameters are important for distance determina¬ 
tion. 







